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<t^ ; Abstract 

^*' . We present a model for the description of orientational relaxation in hydrogen-bonding liquids. 

f-^—^ , The model contains two relaxation parameters which regulate the intensity and efficiency of dissi- 

O ' 

c/5 ' pation, as well as the memory function which is responsible for the short-time relaxation effects. 

^ • It is shown that the librational portion of the orientational relaxation is described by an algebraic 

~ t~^''^ contribution, on top of which more rapid and non-monotonous decays caused by the mem- 

^ I ory effects are superimposed. The long-time behavior of the orientational relaxation is exponential, 

although non-diffusional. It is governed by the rotational energy relaxation. We apply the model 

to interpret recent molecular dynamic simulations and polarization pump-probe experiments on 
> 
^ ; HOD in liquid D2O [C. J. Fecko et al, J. Chem. Phys. 122, 054506 (2005)]. 
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I. INTRODUCTION 

The recent advances in nonlinear ultrafast polarization-sensitive spectroscopy^ 
make it possible to monitor molecular rotation in hydrogen-bonding liquids in real 
tim o^i^i'^i^i^i'''i^i^i-^°i-^^ ■ Due to the enormous complexity of the problem, which is exacer- 
bated by many-body effects and multitudes of the time scales involved, experimental data 
alone are insufficient for understanding the underlying dynamics. Nowadays, molecular 
dynamic has become a standard tool for studying molecular reorientation in liquidsi^. 
Furthermore, one can even use ab initio molecular dynamics (in which the density func- 
tional theory is invoked to describe molecular electronic structure and inter- and in- 
tramolecular forces are calculated on-the-fly)^^^ or centroid molecular dynamics (which 
accounts for quantum effects/corrections)^. On the other hand, there exists a plenty of 
"old" phenomenological models of molecular reorientation in gases and liquids. We men- 
tion the small-angle rotational diffusion mode l^^i^^i^^ , the jump diffusion mode l^^i^^i^^i^^ , 
the friction mode P^i^'^i^^ , the Gaussian cage model^^^, the itinerant oscillator model^, 
along with the more sophisticated memory function approach p8i29i30i3i^32^33 ^^^^ references 
therein), the extended diffusion model o^'^i^^i^^i^'^i^^i^^i"^°i"^^i"^^ , the rotational Fokker-Planck 
equatio n^'^i^^i^^i^^i"^^i"^^i"^^ , the confined rotator model^, the Steele model^, the Keilson-Storer 
model (KSM) ^^i^^i^^i^°i^^i^^i^^ , the fluctuating/stochastic cage mode Pi^'^i^^i^^i^''' , and the gen- 
eralized Langevin equations/normal mode approac h^^i^^i^^i^^ . Furthermore, the model has 
been elaborated^^, which accounts for the effects of rotation- vibration coupling in ensembles 
of hydrogen-bonding molecules on the time-resolved pump-probe signals. Very recently, the 
generalized jump model of water reorientation has been suggested^. The models, of course, 
rely upon a simplified picture of molecular rotation. However, in contrast with molecular 
dynamics simulations, they get a deep insight into physics of molecular reorientation, give 
a clear perception of rotational relaxation and provide us with explicit formulas for the 
pertinent correlation functions (CFs). 

The aim of the present paper is to develop a simple and physically sound model of 
molecular reorientation in hydrogen-bonding liquids. The model is intended to supply ex- 
perimentalists with a simple theory to interpret and to fit their data and to clarify the in- 
terconnection of orientational relaxation and hydrogen bond making/breaking processes. In 
hydrogen-bonding liquids, the angular momentum CFs exhibit pronounced oscillations and 



orientational CFs (OCFs) display a rapid short-time decay followed by a slower (sometimes 
oscillatory) pattern which transforms gradually into a monotonous exponential relaxation. 
By incorporating the proper description of the memory effects into the KSM framework, 
we developed a non-Markovian generalization of thereof, NKSM. Within the NKSM, we 
derived analytical expressions for the angular momentum and energy CFs, as well as simple 
recursive expressions for OCFs. 

The paper is structured as follows. The NKSM is formulated in Sec. II. The explicit 
expressions for the angular momentum CF, rotational energy CF and OCFs are presented 
and discussed in Sec. III. Sec. IV contains illustrative calculations of various NKSM CFs and 
comparisons with the results of molecular dynamic simulations and polarization pump-probe 
experiments on HOD in D2O at a room temperature^. A brief summary of the main findings 
can be found in Sec. V. Appendix A contains the explicit formulas for the calculation of 
spherical and linear rotor OCFs within the NKSM. Analytical expressions for OCFs in a 
particular case of "perfect" librations are obtained and discussed in Appendix B. 

A few words about the notation and conventions, (i) The reduced variables are used 
throughout the article: time, angular momentum and energy are measured in units of 
■\/ 1 / {ksT) , y/IksT and fc^T, respectively. Here ks is the Boltzmann constant, T is the tem- 
perature, and / is a characteristic moment of inertia of the molecule so that r^ = \J I /{kgT) 
is the averaged period of free rotation, (ii) All Laplace-transformed operators are denoted 
by tilde, viz. /(s) = /g°° (itexp{— st}/(t) for V f{t). (iii) Repeated dummy Greek indexes 
imply summation over x, y and z. 

II. THE MODEL 

We start with a formally exact Zwanzig-type master equation, which can be derived from 
the general A^-particle rotation-translational Liouville equation by applying the projection 
operator technique2^>^2iSi 

atp(J,f^,t) = -iA(J,0)p(J,0,t)- / dt'C(J,f2,t-t')p(J,^,0- (1) 

Jo 

Here p(J, fi, t) is the single particle probability density function, J is the angular momentum 
in the molecular frame, O are the Euler angles which specify orientation of the molecular 
frame with respect to the laboratory one. The free-rotor Liouville operator consists of the 



two contributions, 

A(J,f2)=An+Aj, (2) 

which describe, respectively, the angular momentum driven reorientation and the angular 
momentum change during free rotation: 

An = /„ JaLa, Aj = —iSap^Ip JaJpdj^. (3) 

/q, are the main moments of inertia. La are the angular momentum operators in the molecular 
frame. For linear and spherical rotors, Aj = 0. 
The relaxation operator C assumes the form 

C{i,n,t) = C{i)g{t), (4) 

g{t) being the memory function which is normalized to unity, j^ dtg{t) = g{0) = 1. All 
the formulas derived in the present paper are valid for any functional form of g{t). Since a 
simple exponential memory function is known to exaggerate oscillatory effects in the angular 
momentum CF and OCFs (see, e,g. ^25,29^32^50 -^^ ^.j^g two-exponential memory function will be 
adopted for making all specific calculations, viz., 

g{t) = aXi exp{-Ait} + (1 - a)X2 exp{-A2t}, (5) 

s + Ai s + A2 
The parameters Aj regulate the memory effects of the two contributions, and a controls 
their relative significance. If we let both Ai and A2 tend to infinity, then g{t) —>■ 6{t) and 
the Markovian limit is recovered. 

Since molecules are massive inertial particles, the relaxation operator C{3) is assumed 
to be fi-independent^. This is tantamount to the statement that molecular reorientation is 
driven by the (time-dependent) angular momenta, whose relaxation, in turn, is governed by 
operator (^. This assumption is consistent with classical molecular dynamics simulations, 
in which one integrates equations of motion of the kind 

dtD^in) = -lAiiD^in), dtJ = -iAjJ + N, 



D^{Q) being the Wigner D-functionsfii and N being the torque acting on a chosen molecule 
from its neighbours. As is demonstrated below, N is essentially non-Gaussian and non- 
Markovian. 

The operator C can further be represented in the general form^i 

c{j)p{j, n, t) = -z/{p(J, n,t)- f dj'T(j| j')p(J', n, t)}. (7) 

The rate z/ determines the dissipation strength, and the relaxation kernel T obeys the nor- 
malization 



f dJT{3\3') = l (8) 



and the detailed balance 

T(J|J')pb(J')=T(J'|J)pb(J), (9) 

Ps{3) = {27r)-'/'exp{-Jl/{2Q} (10) 

being the equilibrium rotational Boltzmann distribution. 



To proceed further, we adopt the KSM parametrization of the relaxation kerne 



133.64. 



T(J|J')= n ^"(-^^l^a), (11) 

a=x,y,z 

UJalJ'a) = [2vr/a(l - Ta')]"'/' exp{-( J, - ^X)V[^W " ll)]} ■ 

Here the parameters — 1 < 7a < 1 determine the relaxation mechanisms. When 7^ = 1, then 
Ta{Ja\J'a)= ^{Ja — J'a) ^"^^ C'(J) = 0. The Fokkcr-Plauck relaxation operator is recovered in 
the limit 7a— >1, z/— >oo, z/(l — 7a) -^ Va = const. If 7a = 0, intermolecular interactions 
are so strong that they "immediately" restore an equilibrium Boltzmann distribution in 
the molecular ensemble (T(J|J') -^ pb{J))- Therefore, the KSM contains the J-diffusion 
mode l^^i^^i^^i^^i^^i^^i"^°i"^^i"^^ and the rotational Fokker-Planck equation^^iS^iSLiMii^ as special 
cases. 

By letting 7a = —1, one gets Ta(Ja|-'^) = S{Ja + J'a)- Thus the magnitude of the angular 
momentum is preserved but its direction is reversed. This regime (7a ~ — 1, a molecule 
rotates back and forth within the cage formed by its nearest neighbors) is expected to be 
of particularly relevance for hydrogen-bonding liquids. Such physical picture of molecular 



reorientation in liquids is inherent in a number of theoretical approaches, in which the 
influence of the nearest neighbors on the selected molecule is modeled by external potentials 
with several minim a^^i^'''i^'''i^^i^^i^'''i^^i^^i'^°i''''^ , or by fluctuating torques and structuresi^i^L^. 
Within the present approach, the fluctuating cage potential is not introduced explicitly, but 
its influence is taken into account dynamically via kernel (jllj) . 

Before embarking at particular calculations, it is useful to estimate values of the param- 
eters z/, 7a (Eq. (HH)) and Aj (Eq. (jSJ for hydrogen-bonding liquids. Since the molecules 
are assumed to undergo hindered librations, one expects z/ ^ 1 and 7^ ~ —1. The memory 
effects are supposed to be quite significant (Aj ~ 1). 

III. CORRELATION FUNCTIONS 

After the explicit form of the relaxation operator (fTT|) has been determined, the master 
equation ^ can be invoked to calculate any rotational and/or orientational CF of interest. 
To study the evolution of any quantity, which depends solely on the angular momentum, we 
integrate Eq. (Hj) over J7 and obtain the reduced master equation 

dtpiJ, t) = -2Ajp(J, t) - / dt'g{t - t')C^(J)p(J, t')- (12) 
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OCF of the rank j is defined through the Wigner D-functions as follows: 

G\t) =< U' {Sl{t))D\^{^))-^ >= fdJG^{J,t). (13) 

After the insertion of the above definition into Eq. (^ one obtains the following equation: 

/•* 
dtG^iJ,t) = -i{An+Aj)G^iJ,t)- dt'g{t-t')CiJ)G^i3,t'). (14) 

Jo 

Here operator A(J) is determined by Eqs. (J2I) and (jH)), in which the angular momentum 

operators La are replaced by their matrix elements L^ over the D-functions: 

{Li)ki±z{Li)ki = Sk,i^i{U±l)UTl + l)V^\ {L%i = l5ku -J<k,l<j. (15) 

Eqs. fll2j) and (fT^ can be solved numerically in case if a general asymmetric top molecule. 
To make the presentation simpler, we restrict ourselves to the consideration of spherically 
symmetric molecules. The corresponding formulas are much more elucidating and convenient 



to analyse, since the relaxation operator C{3) is described by only two dynamic parameters 
(the intensity, u, and the efficiency, 7, of dissipation) and a spherical molecule possesses a 
single moment of inertia, I. This theory can be applied to asymmetric tops as well. Indeed, 
in the hindered rotation limit {tj ^ 1), a molecule librates back and forth in its cage, and 
every single libration reorients the molecule to a small angle. Using the explicit form of the 
operator A^ (jH)), it is easy to demonstrate that the averaged inertial reorientation angle of 
an asymmetric top around its z-axis equals (t/r^)^ <^ 1, where 



Tr = ./TJik^), r' = {!-' + Iy')/2. (16) 

This corresponds to the rotation of the spherical molecule with the effective moment of 
inertia I. Thus, all the formulas obtained below can be used for asymmetric top molecules by 
making use of the substitution (fTBj) . provided one is interested in the orientational relaxation 
of the tensor with nonzero components along the molecular z-axis. If the quantity under 
study possesses nonzero components along several axes of the main moments of inertia (like 
OH stretch in HOD), the above analysis remains true if the effective moment of inertia / 
is modified accordingly. On the contrary, free inertial rotation of spherical, symmetric and 
asymmetric tops is very differen1>22iSi2i. 

Starting from Eq. (fT^ . it is straightforward to derive the explicit formulas for the Laplace 
transformations of the angular momentum and rotational energy CFs: 

cAt) - ^. CAs) ^ j^, (17) 



'^'=<') = -,4 „■»/ • c-^(^) = —^hr:^- (18) 



;j 'j'w) - (J')' ^,,. L 

Here the rates 

i/j = z/(l-7), UE = iy{l-^^) (19) 

determine the angular momentum and rotational energy integral relaxation times: 

/"OO 

Tj = / dtCjit) = CjiO) = uj\ (20) 
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/■OO 

TE = / dtCEit) = CEiO) = i^e'- (21) 

Jo 



As to the OCFs, the solution of Eq. ()14|) can be given in terms of recursive relationships 
(or, equivalently, in terms of continued fractions) for any value of the relaxation parameters 
V and 7, as well as for any memory function g{t), see Appendix A. 

IV. ILLUSTRATIVE CALCULATIONS 

In this Section, we present and discuss the results of representative calculations of the 
angular momentum CFs (HIj), energy CFs (jIHl), and OCFs ^ - (jSSI) within the NKSM. 



The free rotation period (fTH|) has been taken as r^ = ^/T/Jk^) = 93fs. This corresponds 
to HOD at T = 300 K (4 = 2.63, ly = 1.85, I, = 0.72 a.m.u. x A^). 

We start from the angular momentum ()17p and energy ()18p CFs. If the two-exponential 
memory function (0) is used, the Laplace images Cj{s) and Ce{s) can be inverted into the 
time domain by solving the pertinent cubic equation. The CFs in the time domain are thus 
determined by linear combinations of one real and two complex conjugated exponentials. 
The results of representative calculations are depicted in Fig. 1. A decrease of the second 
memory parameter, A2, causes a typical transformation of the angular momentum CF, which 
reflects a passage from simple to hydrogen-bonding liquids^^. Since CFs Cj{t) attain negative 
values, their actual decay occur at a timescale of several hundreds of femtosecond, which is 
much longer than the integral relaxation time tj = uj^ = 4.65 fs. The rotational energy 
CFs are monotonous and decay at a much longer time scale of 10 4- 15 ps, which is a direct 
manifestation of the librational motion (7 ~ —1, tj <^ te). 

Let us turn to the study of orientational relaxation. To get a qualitative feeling of 
the influence of the relaxation efficiency 7 on molecular reorientation, let us concentrate 
on the long-time behavior of OCFs and consider the Markovian limit, g{t) -^ (5(t). As 
has been established inM^, the KSM predicts that the smaller is 7, the slower is the OCF 
decay. Thus, for a fixed angular momentum relaxation rate z/j, the rotational Fokker-Planck 
equation (7 = 1) predicts the most rapid orientational relaxation, while the limit of perfect 
librations (7 = — 1) corresponds to the slowest orientational relaxation. This is clearly seen 
in Fig. 2. 

As is demonstrated in Appendix B, Eq. (|T4|l can be solved analytically for OCFs in case 



of perfect forward-backward librations (7 = — 1 ) in the hindered rotation hmit (z/ ^ 1): 



k 
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(22) 



On the scale of Fig. 2, the exact solution of Eq. (fT^ in the limit of 7 = —1 and the 
approximate one, which is delivered by Eq. (J22I), are indistinguishable. Thus the in-cage 
librations manifest themselves through a slow algebraic t~^/^ decay of the OCF. This behav- 
ior is caused by the angular momentum reversion (7 = — 1) and has nothing in common with 
long-time hydrodynamic tails of the angular velocity CFs of Brownian particles (see^^iSiZ^ 
and references therein)^^. It is interesting to point out that the OCFs calculated within 
the M-^2iS and E— diffusion models exhibit similar (~ t~^/^ ) long time tails, which are 
commonly regarded as unphysical. Nonetheless, these models were successfully invoked to 
reproduce "experimental" OCFs, which were obtained through the inversion of IR and Ra- 
man spectra into the time domain^iiiSiSiSMi. The present consideration reveals that this 
success is not accidental, since the conservation of the magnitude of the angular momentum 
(M- diffusion) or rotational energy (E- diffusion) mimics the description of in-cage librations 
via the KSM kernel (jllll . which in the limit 7 ~ — 1 conserves both of these quantities. 

Fig. 3 illustrates the influence of the memory effects, which modify the short-time behav- 
ior of the OCFs. The mechanism of this influence is uncovered by the general expressior^ 



G\t) ^ 1 - j(j + 1) / di!{t - t')Cjit'), (23) 

which relates the short-time OCF behavior with the time evolution of the angular momentum 
CF. That is why the parameters which have been used for the calculation of the OCFs in Fig. 
3 and the angular momentum CFs in Fig. la are the same. Generally, the memory effects 
speed up the short-time OCF decay. If the angular momentum CF is highly oscillatory, 
these oscillations show up in OCFs also (compare dotted lines in Figs, la and 3). 

OCFs depicted in Fig. 3 look qualitatively similar to those obtained via computer 
simulationa^ i^i^^i^^i^^i^^ . A more quantitative comparison of the simulated and NKSM OCFs 
is presented in Fig. 4. It depicts the results of molecular dynamic simulations of the first 
(upper dotted line) and second (lower dotted line) rank OCFs for HOD in D2O at a room 
temperature using the SPC/E model for water^ along with the fits obtained within the 
NKSM (full lines). Detail of the molecular dynamics simulation protocol can be found in 
ref.i^ 



The present theory is seen to reproduce the simulated OCFs quite welL We were unable, 
however, to fit the first and second rank OCFs by the same set of the NKSM parameters. 
This is not unexpected: OCFs of different ranks are affected by the cage potentials in 
a different way, which is determined by the potential symmetry^^. The dielectric friction, 
which governs orientational relaxation in polar systems is also known to be rank-dependenli^. 
The values of u, 7, Aj and a which are extracted from the angular momentum CF deliver, in 
fact, certain averaged values. It is nonetheless quite remarkable that the parameters which 
have been used for the calculation of the second rank OCF are identical to those which have 
been used for the computation of the angular momentum CF (Fig. la, full line). This latter 
CF looks very similar to the simulated one^^. 

Fig. 5 shows the normalized anisotropy extracted from the pump-probe signal^ (dotted 
line). It deviates quite significantly from the simulated second-rank OCF (our best fit to this 
OCF is plotted here for the sake of comparison). As has been pointed out i n^i^i^^ , simulations 
normally predict faster, in comparison with experiment, anisotropy decays. This is caused, 
perhaps, by ignoring either the water polarizability^^^ or quantum effecta^^. We are not 
attempting to resolve this controversy here. Note merely that the experimental anisotropy 
can be fitted quite well within the NKSM (full line), with the set of parameters which predict 
more "librational" reorientation (7 is closer to —1) and more pronounced memory effects 
(Aj is smaller). 

Let us return back to Fig. 4. The OCFs calculated within the standard diffusion equation, 

G^(t)=exp{-j(j + l)rjt}, (24) 

which reproduces the long-time limit of the first cumulant formula, 

G^it) = exp{-j(j + 1) / dt'it - t')Cj(t')}, (25) 
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are seen to deviate significantly from the simulated/NKSM OCFs. Despite the simu- 
lated/NKSM OCFs do exhibit the long-time exponential behavior, and despite the rota- 
tional motion is definitely hindered {tj ^ 1), the diffusion equation cannot reproduce these 
OCFs. The failure of the diffusion equation has quite an evident explanation. If we assume 
that the molecules undergo "perfect" librations (7 = —1 within the present theory), then 
the long-time behavior of the OCF in the hindered rotation limit is described by Eq. fl22|) 
rather than by the small-angle diffusion equation (|^1|). Since the actual librations are not 
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"perfect" (7 is close to but less than —1) one expects deviations from Eq. ()22|). As is seen 
from Fig. 6, this is indeed the case. This Figure reproduces the Markovian limits of the 
KSM OCFs from Figs. 4 and 5, along with their counterparts calculated via Eq. (j22|) . The 
short-to-intermediate-time resemblance of the OCFs is quite remarkable. 

On the other hand, the long-time behavior of the simulated/KSM OCFs is seen to be 
exponential and is not reproduced by Eq. ()22|). This hints at a possibility that rotational 
energy relaxation might be responsible for the long time exponential decay of the simu- 
lated/NKSM OCFs. This hypothesis is corroborated by the observation that the fit of the 
simulated and NKSM OCFs of the first and second rank (Fig. 4) delivers different values 
of the relaxation parameters v and 7. However, the quantity v^ = i^{l — 7^) = 0.027 turns 
out to be the same for both j = 1 and 2. This gives the estimated value of 3.4 ps for the ro- 
tational energy relaxation time te- This value correlates with the experimentally measured 
long time anisotropy decay times of 2 -^ 3 p ci^i^i^/i^i^° . 

Note, finally, that a close interrelation between the hydrogen bond dynamics and rota- 
tion dynamics has repeatedly been emphasized in the literature. The NKSM values of the 
"persistence time" of the oscillatory angular momentum CF (Fig. la) and the rotational 
energy relaxation time, ~ 150fs and 3.4ps, correlate quite well with the estimations for the 
continuous and intermitted hydrogen bond lifetimes^i^. This observation is consistent with 
the physical picture of molecular rotation underlying the NKSM approach. The in-cage 
librations, which can cause bond breakings, manifest themselves on the time scale of the 
"persistence time" of the angular momentum CF. On the other hand, the bond rupturings 
are accompanied by subsequent bond reformings, since the molecule does not leave its local 
cage. A "true" cleavage of the bond can occur on the te timescale since, within the NKSM 
description, the rotational energy relaxation occurs due to a hopping to a new position of 
the local equilibrium or restructuring the local potential well. 

V. CONCLUSION 

We have developed a NKSM description of the orientational relaxation in hydrogen- 
bonding liquids. Within the NKSM, molecular rotation is governed by two relaxation pa- 
rameters z/ and 7 (which describe the intensity and mechanism of dissipation), as well as 
by the memory function g{t) (jSj), which is responsible for the short-time dynamics. Alter- 
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natively, the relaxation parameters v and 7 are uniquely determined through the angular 
momentum and energy relaxation times tj ()20|) and te ()21|) . Once a set of the parameters 
is selected, the NKSM allows to calculate any rotational CF or OCF of interest. Keeping in 
mind a considerable success of the KSM in reproducing molecular reorientation in gases^^^, 
the results of the present work demonstrate that the NKSM can be used for the description 
and interpretation of the orientational relaxation in a condensed phase, from rarefied gases 
with binary collisions, through dense fluids to hydrogen-bonding liquids. 

The NKSM suggests the short-time relaxation of the OCFs in the hydrogen-bonding 
liquids is described by an algebraic ~ t~'^' ^ contribution. This algebraic behavior is modified 
by more rapid and non-monotonous dynamics, which is induced by the memory effects. The 
long-time decay of the OCFs is exponential, although non-diffusional. It is governed by 
the rotational energy relaxation time, te- Our results are contrary to standard belief that 
the angular momentum CF determines molecular reorientation in the hindered rotation 
limit, and the first cumulant expression, Eq. ()25|1 . delivers the leading contribution into the 
OCF. Our results indicate, that knowing Cjif) is not enough to predict OCFs for hydrogen- 
bonding liquids, since the long-time behavior of OCFs is governed by the rotational energy 
CF, Cs(t). 

It is a conventional practice to fit various experimental or simulated CFs via a linear 
combination of several real (if the CF decays monotonously) or complex (if the CF exhibits 
oscillatory behavior) exponents. According to the present analyses, the OCF in hydrogen- 
bonding liquids contains an algebraic ~ t^^/^ contribution. This finding suggests that the 
following fitting formulas for the angular momentum CF, 

Cj{f) = fli exp{— z/it} + exp{— z/2t}(a2 cos(f2t) + 03 sin(r2t)), (26) 

and for the OCF, 

G^{t) = hi exp{-z/it} + exp{-z/2t}(62 cos{VLt) + 63 sin(fit)) 

+64 exp{-z/3t} [1 + h^tY^'"^ + 66exp{-z/4t}, (27) 

can be more physically motivated (a,, 6j, z/j and Vt being certain real- valued parameters). 
Indeed, Eq. ()26|) is nothing else than a formal solution of Eq. ()17p in the case of two- 
exponential memory function (jSJ. As to the OCF (jTfj) . the first three terms with coefficients 
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bi, 62, &3 describe the angular momentum induced short-time rapid decay and oscillations. 
The term which is proportional to 64 is responsible for the algebraic contribution, and the 
last term governs the long-time exponential decay. 
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VI. APPENDIX A. RECURSIVE RELATIONS FOR LINEAR AND SPHERICAL 
TOP OCFS 

After being Laplace transformed, Eq. (jl4|) reads: 

-Pb{J) + s&{3,s) = -tiAn + Aj)G'iJ,t)-ug{s){G'{J,s)- f d3'T{J\J')G'{J', s)}. (28) 

Since Eq. (|28|) depends on s parametrically, the method of its solution in the Markovian 
limit {g{s) = 1), which has been been developed i n^5^46^5i^52 ^ jg directly applicable to the 
present case also. One has merely consider the complex quantity lygis) as the generalized 
relaxation rate. We therefore present the final expressions for the calculation of the first and 
second rank OCFs. 

Their Laplace images of the spherical top OCFs can be calculated via the formula 

G\s) = {l + 2bo)/s. (29) 

For j = 1, the coefficient 60 can be retrieved from the simple three-term recursive formul a^^i^^ 
1 _4m + 10 (2m + 3 2m + 2 1 m 

-OmO — Om+1 — i 1 \- Cm ( Om -\ Om-1, (30j 

■5 0"m+l t ^m+1 (^m ) (^m 

6mQ being the Kronecker delta. The value of Bq for the second rank OCF can be extracted 
from the system of coupled recursive relations for the coefficients bm and drrr^: 

1^ _ r 6 1 12 

-OmO — — S \- Qm ( bm -\ Om+1 (jlj 



4m f8m-14 8m - 31 1 -16m^ + 74m + 12 
dm-i + < \ } dm H 7 — T c 

CTm [ CTm 0"^+! j [m + Ipm+l 
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m f 2m — 1 2m — 1 1 , 4m , ^ 

= — bm-i - < \ > b-m H bm+i (32) 

5m J 14 4 + 10m . 1 , 4m(9 + 5m) 

H d-m-l ~ S 1 I" Cm r dm + 7 ""TT "m+1- 

0"m lO-m 0"m+l J (m + 1)0-^+1 

Here 

a^ = s + i^~g{s)i^ - I'n, Cm = s + i^~g{s)i^ - 7'"^+')- (33) 

The solution of the recursive relations for the first rank OCF ^^ can be expressed in the 
continued fraction fornij^^, while the solution of Eqs. (jHlll and ()H2|1 for the second-rank OCF 
can be given in terms of the matrix 2x2 continued fractions^. 

The first and second rank OCFs for linear rotors can be evaluated very similarly, through 
the simple three-term recursive formulas. Namely, the Laplace images of the OCFs can 
also be computed via Eq. (|^ . For j = 1, the coefficient bo must be determined by the 
formul a?^'^^ 

1 4m + 8 r2m + 2 2m + 2 ^1, m 

-OmO = bm+1 — i 1 \- Cm ( bm -\ Om-1, 

while the second-rank OCF {j = 2) can be computed through^ 

3 _ 16m + 32 rSm + lO 8m + 6 1 4m 

-OmO — bm+1 — S 1 \- Qm } bm -\ Om-1- 

S CTm+l L ^m+1 <^m J CTm 

VII. APPENDIX B. ORIENTATIONAL RELAXATION IN CASE OF PERFECT 
ANGULAR MOMENTUM REORIENTATION 

If we neglect the memory effects {g(t) — > S(t)) and put 7 = — 1, then Eq. (fT^ can be 
solved analytically: 



G^(i) = ^— Yy^ dJpBiJ) (34) 



1 + exp{— i/t} y^ [(1 + u/uJk) expjcjfct} + (1 — i^/i^k) exp{— cut}] 



3 

\ fc=i 

The frequencies are explicitly defined as follows 



tffc 



Vz/2 - A;2J2. (35) 
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Several important properties of Eq. (jMj) are to be discussed. OCF (|Mj) possesses a stationary- 
asymptote: G^{t —>■ oo) = (2j + 1)^^, which is identical to the free OCF asymptote. This 
can be easily understood: reversion of the angular momentum is equivalent to the reversion 
of the sense of molecular rotation. Therefore a sequence of forward-backward rotations is 
equivalent to a single free rotation. This asymptote is solely caused by dynamic effects 
(in-cage librations, compare with®) rather than by external potentials (see, e.g.^^^^^). The 
librational motion itself is caused, of course, by in-cage potentials but they do not enter 
explicitly into our analysis. If z/ ^ 0, Eq. ()34|) reproduces the free spherical top OCF. If we 
take the opposite limit u -^ oo, then G^{t) -^ 1 since a large number of small-angle forward- 
backward rotations causes no net reorientation. In the hindered rotation limit {tj ^ 1) Eq. 
(IMj) reduces to (J22D. 

Using the explicit form of the free linear rotor OCF— one can easily derive the linear 
rotor counterpart of Eq. 



l + —t 



2 1-1 

(36) 



fc=i 

^kmif^) being the reduced Wigner function^. It is well known that the exponent d of the 
long-time hydrodynamic tails t"*^/^ of the angular velocity CFs of Brownian particles is 
determined by the dimensionality of the rotation space: d = 3 for any spherical, symmetric 
or asymmetric top while rf = 2 for a linear rotor (2^i2^ and references therein) . This general 
statement holds true in the present case also, and OCF (jHHjl possesses a ~ t~^ tail. This 
means that the algebraic contribution to OCF ()36|) decays slower than its t^'^^'^ counterpart 
in the spherical top OCF (j221)- Furthermore, (ioo(^/2) = for odd j. Thus, the odd-ranked 
OCFs do not possess the stationary contribution and decay faster then the even-ranked 
OCFs. 
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FIG. 1: Angular momentum (a) and energy (b) CFs for z^ = 16, 7 = —0.99915 (that is vj = 32 
and ve = 0.027), a = 0.2 and Ai = 0.7. The dashed, full and dotted lines correspond to A2 = 1000, 
10 and 3, respectively. On the scale of the figure, the rotational energy CFs for A2 = 10 and 3 are 
indistinguishable. 
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FIG. 2: The second rank OCFs in the Markovian hmit for the relaxation rate rate uj = 32. From 
bottom to top, the curves correspond to 7 = 1, —0.99, —0.999, —0.9999 and 1. 
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FIG. 3: The influence of the memory effects on the second rank OCFs. I'j = 32, 7 = —0.999, 
a = 0.2 and Ai = 0.7; A2 = 1000 (dashed lines), 10 (full lines) 3 (dotted lines). The dash-dotted 
curve depicts the OCF in the Markovian limit. 
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FIG. 4: Comparison of the simulated OCFs with those calculated within the NKSM. The dot- 
ted lines correspond to the first (upper curve) and second (lower curve) rank OCFs which were 
simulated for HOD in liquid D2O at a room temperature^. The black curves are computed for 
uj = 46, 7 = -0.999415, a = 0.2, Ai = 0.4, A2 = 10 (j = 1); vj = 32, 7 = -0.99915, a = 0.2, 
Ai = 0.7, A2 = 10 (j = 2). The OCFs calculated via the diffusion equation (pi)) are depicted by 
dashed lines. 
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FIG. 5: Comparison of the experimental OCFs with those calculated within the NKSM. The 
dotted line reproduces the best fit to the experimental anisotropy decay^. The solid black curve 
shows the second-rank OCF computed within the NKSM for uj = 63, 7 = —0.99995, a = 0.2, 
Ai = 0.3, A2 = 10. The dashed line reproduces the second rank NKSM OCFs from Fig. 4. 
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FIG. 6: Elucidation of the algebraic contributions into OCFs. For each couple of the curves, 
the upper one is calculated via Eq. (|22|) and the lower one is computed within the Markovian 
limit of the NKSM. The dotted (j = 2, vj = 32, 7 = -0.99915) and dashed {j = I, vj = 46, 
7 = —0.999415) curves correspond to the best-fit NKSM OCFs from Fig. 4, and the solid curves 
(j = 2,uj = 63, 7 = -0.99995) correspond to the best-fit NKSM OCFs from Fig. 5. 
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